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Boson stars in zero-, one-, and two-node equilibrium states are modeled numerically within 
the framework of Scalar- Tensor Gravity. The complex scalar field is taken to be both massive 
and self-interacting. Configurations are formed in the case of a linear gravitational scalar coupling 
(the Brans-Dicke case) and a quadratic couphng which has been used previously in a cosmological 
context. The coupling parameters and asymptotic value for the gravitational scalar field are chosen 
so that the known observational constraints on Scalar- Tensor Gravity are satisfied. It is found 
that the constraints are so restrictive that the field equations of General Relativity and Scalar- 
Tensor Gravity yield virtually identical solutions. We then use catastrophe theory to determine the 
dynamically stable configurations. It is found that the maximum mass allowed for a stable state in 
Scalar- Tensor Gravity in the present cosmological era is essentially unchanged from that of General 
Relativity. We also construct boson star configurations appropriate to earfier cosmological eras and 
find that the maximum mass for stable states is smaller than that predicted by General Relativity, 
and the more so for earlier eras. However, our results also show that if the cosmological era is early 
enough then only states with positive binding energy can be constructed. 
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I. INTRODUCTION 



General Relativity describes very well the gravitational interaction. However, it still has competitors, one of which 
is Scalar-Tensor Gravity (see Will Q for a review). Besides the metric, Scalar- Tensor Gravity has an additional scalar 
field. Motivation for this additional field can be found in Dicke's |^ discussions about arbitrariness in the measurement 
of lengths, times, and masses. There is also the fact that dilaton gravity is the low string tension limit of string theory 
(see, for instance, Damour and Polyakovjsl). Motivation for the theory itself, is the attractor mechanism discussed 
by Damour and Nordtvedt j|] (see also [^H])- They demonstrate that the homogeneous and isotropic expansion of 
the universe forces a series of Scalar- Tensor cosmological solutions to be virtually indistinguishable from General 
Relativity solutions at late cosmological times. This implies that today the weak-field differences between Scalar- 
Tensor Gravity and General Relativity are small. It is the purpose here to see if significant strong-field differences 
between the two can be generated for equilibrium state boson stars. 

A boson star consists of massive scalar particles (i.e., excitations of a complex scalar field) that are held together by 
gravity, a self-interaction, or both. They were first discussed by Kaup ||^ and then by Ruffini and Bonazzola ||^ (for 
a thorough review, see Liddle and Madsen [^). They can form stable configurations having negative binding energy. 
If we include a self-interaction, then their maximum allowed stable mass comes to the order of a solar mass [|lO|. It 
is even speculated that they are a form of dark matter and that they can be created during a phase transition in the 
early universe (see Frieman et al pT|]). Unlike boson stars formed today, those formed in the early universe would 
have done so in an era when even the weak-field differences between Scalar- Tensor Gravity and General Relativity 
are significant. 

In terms of their equilibrium characteristics, there are infinitely many different boson star configurations for a 
given central density, depending on how many nodes, or zeroes, there are in the matter scalar field. Furthermore, 
the ground state, or zero-node configuration, has the lowest ADM mass, the higher-node configurations then have 
successfully higher masses. With this in mind, there are three different ways that we try to generate non-trivial 
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Scalar- Tensor effects: (i) build zero- and higher-node configurations for three different values of the matter scalar 
field self-interaction, (ii) build sequences of zero-node solutions, where different sequences have different values for the 
self-interaction, and (iii) repeat (i) and (ii) for two different Scalar- Tensor theories. 

The two varieties of Scalar- Tensor Gravity considered here use the linear Brans-Dicke coupling and the quadratic 
coupling used by Damour and Nordtvedt Q . Experiments performed in our solar system have placed a lower bound on 
the slope of the Brans-Dicke linear coupling . We will only consider Brans-Dicke boson stars that are consistent with 
this lower bound. As for the quadratic coupling, we impose boundary conditions such that the gravitational scalar 
can have an asymptotic value equal to the ambient, cosmological value existing at the present time, as well as earlier 
times. The ambient value is derived from a cosmological solution that exhibits the attractor mechanism of Damour 
and Nordtvedt [Q. When the present-day ambient value is used, then the present-day observational constraints are 
satisfied. 

Boson stars in Scalar- Tensor Gravity have been investigated previously by Gundersen and Jensen JTst , who used 
the Brans-Dicke coupling, and by Torres Q, who looked at three different classes of couplings. Both discussions 
considered only zero-node configurations. Here we will consider solutions containing up to two nodes. Furthermore, 
neither of the previous discussions considered the question of stability. Here, we will apply catastrophe theory to the 
binding energy to determine the change of stability for the zero-node solutions. Finally, we, like Torres, will look at 
boson stars formed in earlier eras of the history of the universe. 

We will show below that stable, zero-node boson star solutions in Scalar- Tensor Gravity are, for all practical 
purposes, indistinguishable from those of General Relativity in the present cosmological era. We also show that the 
same holds for the higher-node solutions, regardless of the strength of the self-interaction. It is for this reason that 
we construct zero-node boson star configurations appropriate to earlier cosmological eras, so as to determine if the 
indistinquishability of the solutions is a generic feature of boson stars or a consequence of the cosmological conditions. 
We shall see that it is the cosmology that is predominate, in agreement with the suggestion of Torres [|4| that masses 
can decrease from their General Relativity values for earlier cosmological eras. But moreover, we use catastrophe 
theory to show that the maximum, stable mass decreases as the cosmological time decreases. However, when the 
cosmological era is early enough, we find that only positive binding energies result so that no stable configuration can 
be constructed. 

The fact that negative binding energy states can exist at all is a priori a non-trivial result. Such states can be 
realized within General Relativity because a boson star's typical mass is M ~ l/Gm and its total number of 
particles is Np ~ 1/Giv? — where G is the Newton gravitational coupling and m is the individual mass of the scalar 
particles making up the star — and so both are of comparable value in the binding energy M — mNp. However, one 
way in which Scalar- Tensor Gravity differs from General Relativity is in how it locally varies the Newton gravitational 
coupling (see, for instance, Bekenstein jlSj and Bekenstein and Meisels and references therein). The inherent 
non-linear features of the Scalar- Tensor field equations suggest that a spatially varying G, which at a very basic level 
controls the strength of the gravitational interaction between two particles, should lead to qualitative differences in 
binding energy between General Relativity and Scalar- Tensor Gravity. We shall derive below a necessary condition 
for negative binding energies to exist. 

In the next section we will introduce the field equations for the gravitational fields and the matter. In Sec. IH we 
write down the equations in spherically symmetric form. Similar to the approach of Friedberg et al |p^ , we write 
down a total energy functional from which the field equations can be derived. It is also here that we discuss the 



boundary conditions and the numerical technique. In Sec. IV we give a detailed account of the solutions containing 
up to two nodes for a typical central value of the matter scalar field. In Sec. we investigate the stability of the 
ground state configurations throug h the use of catastrophe theory. In Sec. ^ we discuss boson stars formed in earlier 



cosmological eras. Finally, in Sec. VII we make some concluding remarks. Units such that c — 1 and h = 1 will be 
used. This implies that the scalar field mass m is an inverse length (actually, the inverse Compton wavelength of the 
scalar particles) and the bare gravitational constant G* has units of length squared. 

II. THE FIELD EQUATIONS 

The action for our system of Scalar- Tensor Gravity coupled to a self-interacting, complex scalar field in the physical, 
"Jordan-frame" is 




B.T. , (2.1) 
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where B.T. represents the boundary terms that can be added to subtract out the second-order derivatives coming 
from R 1^,0 (the expUcit form being identical to that used by Friedberg et al ||l^] and wiU be given below in the 
"Einstein-frame"). The gravitational scalar is (j) and uj{4)) is the "Jordan-frame" coupling of to the matter. The 
complex scalar (with its complex conjugate being t/j'^) has mass m and is self- interacting through the potential 

There is an alternative representation of the action above, for the so-called "Einstein-frame." The transition to this 
frame is effected by the conformal transformation 

.9m- = , (2.2) 

where 

= G,e2"('^) (2.3) 

and a{(p) is the functional transformation from to the "Einstein-frame" gravitational scalar ip. The relationship 
between a;(0) and a(ip) is obtained from 



where 

The action in the "Einstein-frame" is thus 
1 

IBttG 



= (2w + 3)-^ , (2.4) 



a(^) . - . (2.5) 



B.T. . (2.6) 



It does not deliver exactly General Relativity because the metric g^^ is not the true, physical metric that encodes the 
distance between spacetime points. However, the "Einstein-frame" does deliver equations that are similar enough to 
General Relativity that we will use it for our calculations. 
The "Einstein-frame" stress-energy tensor is 



T — 



The gravitational field equations for g^^, and ip are 

G^^. = SvrG,!),^ -f- 2d^ipd,,ip - drfd''fg^,^ (2.8) 

and 

V^V'"(p = -A-KuT , (2.9) 
where T is the trace of the stress-energy tensor. The matter field equations are 

V,V"?At + 2aa,^^9V = e'"^^^ (^"^V^ + 2^^^ (2.10) 

and 

V^Vi) + 2adr'ipd^(p = e^''^'^'> (^m^ip + 2^^ • (2.11) 

The coupling function a(ip) is a priori unknown. There are, however, some theoretical reasons to motivate explicit 
forms. Furthermore, o nce a parti cular form for a{(p) is taken, there are experimental constraints p],p2| that can be 
imposed (using Eqs. ( |3.18[ ) and ( 3.19| ) below). The two forms for a{ip) that will be used here are the Brans-Dicke 
coupling 
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a{^) = J" (2.12) 
\J2lobd + 3 

and the quadratic coupling 

aCv') - ^ VL) , (2.13) 

which was the particular form considered by Damour and Nordtvedt (except for the additive constant). The term 
LPao represents the asymptotic value of the gravitational scalar field. It is known from Solar System observations that 
> 500; likewise, it is known from observations of binary pulsars that k > —5 |2^ . 
As for the complex scalar field, we will be considering a self-interaction of the form 

y(^t^) ^|(^t^)\ (2.14) 

The strength of the self-interaction. A, will be taken to be positive. Because the potential V is a functional of 

then it preserves the global U(l) gauge symmetry (■(/; e'^°'ip, where ct is a constant) present in the theory. This 
symmetry results in a conserved current, whose explicit form in the "Jordan- frame" is 

> = le-^'^f^^^'^" {iPd.ijj^ - i^^d,^) . (2.15) 
This conserved current leads to a conserved charge, which is Np^ the number of particles making up the star (cf., Eq. 



(3.11) below). 



III. THE EQUILIBRIUM STATE EQUATIONS 

A. The equilibrium state field equations 

The spacetimes considered here are spherically symmetric and static, with the "Einstein-frame" metric taking the 
form 

ds^ = -N^{r)df + A{r)dr'^ + [dd^ + sin^edcj)^] . (3.1) 
The gravitational scalar, which is real, is assumed also to be spherically symmetric and static: 

^ = ^{r) . (3.2) 
As for the matter scalar field, Friedberg et al p9| show that the minimum energy configurations arc those for which 



-iilt 



$(r) , 



(3.3) 



where is real and positive. Their proof (see the App endi x in ]l9[) also goes through for scalar-tensor gravity (making 
use of an energy functional defined as below in Eq. ( ^.1C| )), and so we will ta ke ip to have this form. 

To write down the field equations for this system we could take Eqs. (3.1-^^3) and insert them into the full field 
equations written earlier. However, it is more useful to insert Eqs. (3.1-3.3) directly into Eq. ( ^.(j ), and then derive 
the field equations from the reduced action. Let Lmatter represent the reduced "matter" Lagrangian (which includes 



the contribution due to ip) , Lgrav the reduced gravitational Lagrangian, and L = L^, 



i.e., S = J dtL where S is the total reduced action. Then for our particular system we find 

Lmatter = ' dr N ^/A{u + v-w) 
Jo 

and 

Lgrav = 



Lgrav the total Lagrangian, 





n ^'1 


1 


TV' 




1 dr N - 2 


1 + r — 




1 + 2r — 


) 






^71 


N 





(3.4) 



(3.5) 



where 
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w 



(3.6) 



U = 2^G, ( 1 + -^$2 ] e4a(^)<j,2 



(3.7) 



and 



(3.8) 



Note that the gravitational Lagrangian does not foUow just from the Hilbert action, it also includes the contributions 
due to the boundary term B.T. of Eq. (2.6): 



2Nr^ f N' 2 



A \N r 



ANr 



(3.9) 



r=0 



We can obtain the total energy functional E (i.e., Hamiltonian) from the Lagrangian via the standard canonical 
transformation that replaces ip with its conjugate momentum: 



E = QNp ~ L 



m 
G 



dr r^NVA {w + u + v) - Lg 



* Jo 



where Np is the conserved total particle number (in the physical "Jordan- frame" ) and is given by 



2m2 , ^ 



(3.10) 



(3.11) 



By construction it is the case that 



dE 



_ dM 



(3.12) 



where M is the "on-shell" total mass-energy (given below in Eq. ( 3.2C )), and Np is also "on-shell." 

There are now two ways to derive the field equations: one can vary L keeping Q fixed to get them, or one can vary 
E keeping Np fixed. Either way gives the same results, which are 



A' 1-A 



+ 2m^rA {w + u + v) 



(3.13) 



= h 2m rA iw — u + v) 

N r 



(3.14) 



r'^N(p' 



2e 



2a{ip) 



n 



2 1 



i^N^ 



^ ['J'f + A^'"^"^1>'') - (3-15) 



Am? 



and 



m'^N-^ 



A 



(3.16) 



5 



B. Boundary conditions 



The boundary conditions for this system of equations must take into account three things: (i) the solutions must 
be geometrically regular at the origin; (ii) the solutions must yield an asymptotically flat spacetime; and (iii) the 
solutions must take into account the cosmological input for both the coupling a{Lp) as well as Lp. 

Geometrical regularity at the origin means there is no conical singularity, i.e., the proper radius divided by the 
proper circumference should reduce to 2tt at r — 0. This implies that A{0) = 1. Also, to maintain regularity in the 
field equations as r ^ 0, we impose that $'(0) = and (p'{0) — 0. 

For a purely technical reason (to be discussed below), we desire solutions that are asymptotically flat in both the 
Jordan and Einstein frames. That is, we want both g^i, and to reduce to the flat spacetime metric at spatial 
infinity. The implication of this is that the value of ipoc = y^{oo) must be such that a(ipac) = 0. The other outcome is 
that both N and A approach one at spatial infinity. 

For the Brans-Dicke coupling, <I>c = <&(0) and (poo are the only freely specified field values. The value of iV(0) is not 
specified freely, but rather is determined so that $00 = ^"(00) = 0. The value of (p at the origin is not specified freely; 
it must be determined in such a way that the solution for tp goes to ipca at spatial infinity. We will use the freedom 
to add an arbitrary constant to the Brans-Dicke coupling a{ip) (to be discussed in more detail below) so that all the 
solutions we consider have ip^o — 0. 

For the quadratic coupling, again <i>c can be specified freely, but we use the cosmological model of Damour and 



Nordtvedt Q to calculate 



o^R f, 3 \ . ( 5 ISk , ^ /8k 



K \ 8k \4V3 V3 



1 - - e-3p/4sin U - lp + arctanW^ - 1 , (3.17) 



where an ~ 1 is the value of a for the universe at the end of the radiation dominated era and p is a measure of the 
time since this era (i.e., today, p ^ 10); in terms of the cosmological redshift z we have z ~ 6^""^ — 1. Note that we 
must take k = da/d(p > 3/8. 



Damour and Nordtvedt used a form of the scalar-tensor coupling without an additive constant (cf. Eq. ( 2.13 )). 
Fortunately, it can be shown that putting in such a constant does not change the solution above for (poo ■ Furthermore, 
the addition of such a constant will not affect the predictions of Damour and Nordvedt for the PPN (Parametrized 
Post-Newtonian) parameters 7 — 1 and /3 — 1 or the time rate-of-change of the Newton coupling, since each of these 
only depend on the present day values for a and k: 



and 



/?-l--^«(7 + l)(7-l) • (3.19) 



Observational constraints (see Will ||l|, and references therein) yield I7 — 1| < 2 x 10 ^ and |/3 — 1| < 2 x 10 



C. ADM Mass 



A consequence of determining (p at spatial infinity such that 0(9300) = is that the ADM mass M is given by 

/>OC 

G,M = mM drr^{w + u + v). (3.20) 
Jo 

This can be shown as follows: the "Jordan-frame" ADM mass Mj is given by 

G,Mj =. lim J (1 - l/~grr) . (3.21) 

r — ^00 Z 

The similar "Einstein-frame" ADM mass Me is 

G,Me = lim J (1 - l/grr) ■ (3.22) 

r — >oo Z 

However, since grr = ^^^'^grr and a{ipoo) = 0, then the limits on the right-hand-sides are equal and therefore 
Mj = Me = M. Eq. (|^) follows using the integrated form of Eq. ( ^.13| ). Furthermore, it can be shown that the 
value for M one obtains from Eq. ( |3.2C| ) is the same as that delivered by the "on-shell" value for the total energy 
functional E of Eq. ( |3T0| ). 
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D. Necessary condition for negative binding energy states 



If the condition 



2 > n/m 



(3.23) 



is sati sfied , tlien nega tive binding energy states can be constructed. This can be estabhshed by first adding together 
Eqs. (3.13) and (3.14), which results in 



(in [iV^A])' = Am^rA{v + w) 



(3.24) 



The right-hand-side of this equation is positive-definite so that (In [A^^^])' > for all r. Our choice of coordinates 
and boundary conditions dictate that 



hm In [N^A) = 



(3.25) 



Therefore, In (A^^^) < for all finite r and thus N'^A < 1. Because the integrand of M is positive definite (for 
positive A), and N^/A < 1 in the integrand for Np, then 



mNp < 2{m/n)M 



(3.26) 



Now, a negative binding energy state is defined by the condi tion M — mNp < 0. Using the inequality in Eq. ( 3.26| ) 
we see M < mNp < 2{m/il)M , and hence the result in Eq. ( ^.23 ) follows. 

Actually, the preceeding discussion does not give the tightest limit on the ratio fl/m. An even tighter constraint 
on Q/m can be found by looking at the asymptotic form for $ as r gets very large. In this limit, the field equation 
for $ becomes (assuming also that is negligible compared to <&) 



$" - (m^ - r2^) $ w , 

the solution of which is $ ^ cxp (— — Jl^ r). Thus, the stronger constraint is that fi/m < 1. 



(3.27) 



E. Rescaled equations 



We will take advantage of scale-invariances of the field equations to redefine some of the fields, parameters, and the 
radial and time coordinates: 



X = mr , v/47rG,$ -> $ , mN/n -> N , A/AnG^m'^ A , nt/m -> t 
The field equations become (' = d/dx) 

A' 1~A 



A X 



+ 2xA {w + u + v) 



(3.28) 



(3.29) 



2N' 1 - A 



N 



2xA {w — u + v) 



(3.30) 



(3.31) 



2.2jyg2a(v>)^/\ ' 



where now 



2N^ ' 



(3.32) 



(3.33) 
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2 
and 

v = ^{e'''^''H<^>f + yf) ■ (3.35) 

The rescaling does not change the asymptotic value of A — it stiU becomes one at spatial infinity — but it does change 
that of N, which is now 

lim N{x) ^ m/n . (3.36) 

X — *oo 

Thus boson stars with negative binding energies will have asymptotic values for (rescaled) N that will never be lower 
than 1. The total mass and total particle number are also changed to 

A/ = / dxx^{w + u + v) ^ (3.37) 

and 

2 P°° 1 
N.^^l dxx^NVAu.^^. (3.38) 

The integrands, as well as the integrals themselves, are dimensionless. Hence, it is the factors in front that determine 
the typical values for M and Np, and they are those discussed in the Introduction. 

There is one more rescaling that can be done that has no analog in General Relativity, and that is an invariance 
of the field equations if an arbitrary constant is added to the scalar-tensor coupling. If we simultaneously do the 
rescaling 

e'^x-^ X , e=$ $ , e^N N , e^A^ A (3.39) 



on the variables defined by Eq. ( 3.2S ) and let a{(p) + c ^ then the field equations remain unchanged. It is for 

this reason that we can maintain all generality and still have boundary condtions for the universe today such that 
a(Voo) = 0. 



F. Numerical technique 



We have extended the code originally developed by Seidel and Suen [g3j for General Relativity. It is based on a 
fourth-order Runge-Kutta algorithm and now solves the Scalar- Tensor Gravity differential equations ( 3.29| )-(3.32). As 
mentioned earlier, our system requires a two parameter search to find a solution that satisfies the boundary conditions 
for $oo and ipoo that were given earlier in Sec. Ill C . Operationally, we choose a central value of the scalar field <I>c 



first together with a guessed central value of the gravitational scalar field (p{0), and integrate out to large radii for 
different values of N{0). We then check if the resulting ipoo is close to our expected boundary value. In order to judge 
the convergence of the matter scalar field $oo, we set the tolerance to 10^^", which means an asymptotic value for <& 
is convergent if it is less than this tolerance. 

As readers will find lat er in F igs. 1-4, the gravitational scalar field falls off more slowly than the matter scalar 
field. The field equations ( 3.29| )-( 3.32 ) imply that the asymptotic behavior is (p ^ ipoo + C/x, where C is a constant. 
Therefore, at the numerical boundary, say x = Xend, we set the expected boundary values for Lp{xend) as 



ip(Xe 



Po 



c 



where 



r--x^ ^ 



dx 



(3.40) 



(3.41) 



If the computed ip{xend) is not the expected value, then we change <^(0) and repeat the whole procedure. We set the 
tolerance to judge convergence in ipoo as 5 x 10~^. 
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The key distinctions between the matter scalar field with no nodes, and one with nodes, is the zero- node field only 
has an extremum at the origin, and approaches zero asymptotically, whereas all higher-node solutions have as many 
new extrema as zeroes in the matter scalar field. Hence, the algorithm that constructs higher-node fields searches 
not only for zeroes, but also for local extrema. Convergence towards a solution satisfying the boundary conditions 
is effected by applying the same two tolerances given above. In particular, when we search for a solution having 
ri-nodes, then we impose the tolerance after the matter scalar field has passed through the (n + l)-th local minimum 
or maximum. 

We have checked our code in three different ways: (i) Produced solutions in the large coupling limit and verified 
that they are identical with original General Relativity solutions; (ii) Replaced a{ip) with a{ip) + c and confirmed, for 
a given set of boundary condtions, that we obtain the same sequence of results; and (iii) Calculated the total mass 



M three different ways (using eqs. (3.2C-3.22)) and obtained consistent results 



IV. THE EQUILIBRIUM CONFIGURATIONS 

The main goal here is to see if significant Scalar- Tensor Gravity effects can be generated in equilibrium state 
boson stars. We strive for this goal in three ways: First, specific equilibrium configurations are constructed for a 
representative choice of the central value of the matter scalar field. We produce different configurations by increasing 
the node-number in the matter scalar field, increasing the strength of the matter scalar field self-interaction, or 
both. Second, we produce three sequences of zero-node configurations, where an individual member of a sequence is 
specified by the central value of the matter scalar field. The sequences themselves are distinguished by the strength 
of the matter scalar field self-inter action. Third, we repeat the previous two steps for the two different Scalar- Tensor 



couplings given in Eqs. ( |2.12| - p.l3D as well as for General Relativity. 

We have in Figs. 1-3 typical plots of the radial dependence of the matter and gravitational fields for the Brans-Dicke 
coupling. Figs, la, lb, and Ic contain plots of $ versus x for n = 0, 1, 2, where n represents the number of nodes in 
the matter scalar field, and A = 0, 10, 100. The corresponding plots for TV and A are given in Figs. 2a, 2b, and 2c and 
those for (p in Figs. 3a, 3b, and 3c. All the configurations are for $c — 0.15, which will b e seen in the next section to 



correspond to stable zero-node states, and ljbd = 600, which can be seen from Eq. ( 3. IS ) to be consistent with Solar 
System constraints. The same set of plots for boson stars in General Relativity are essentially indistinguishable from 
those presented here. 

Notice that the effect of the nodes in the matter scalar field is to slightly ffatten in a small region around the 
X values where the nodes occur. For example, in Fig. 2c, we can see that A has two local extrema produced for 
each node in the scalar field. The local minimum occurs at precisely the value for x where the node occurs. The 
gravitational scalar (p behaves similarly to N near the nodes. This behaviour for each of the fields is a reffection of 
what is occuring in the energy density: it becomes nearly zero in the small region around each node. Thus, instead 
of the steady accumulation in mass that usually occurs as one moves radially outward from the center of a star, we 
have near each node essentially no accumulation. The metric function A(x), for instance, then behaves like its black 
hole counterpart and decreases as a node is approached (since the mass enclosed within x remains nearly constant 
near each node). After passing through a node, the mass will start to accumulate again so A{x) will start to grow. 
Similar remarks apply to the other fields. 

Initially, one might have expected the gravitational scalar field to have had some effect near the nodes. After all, 
the mass enclosed inside a given rad ius d ep ends on the quantities wlx), u{x), and v{x). They, in turn, depend on the 



gravitational scalar field (cf., Eqs. ( j3.33 )-( p.35| )). It is because the gravitational scalar remains essentially constant 



(and that the realistic boundary conditions coming from the observational constraints force it to be small) that it has 
no significant effect on the accumulation of mass near the nodes. This, in turn, is due to the fact that the matter 



scalar field serves as a direct source for the gravitational scalar field (cf., Eq. (3.31)). 

Using Tables I and H we can compare the masses, particle numbers, radii (which is defined to be the radial coordinate 
value from the star's center containing 95% of the mass), and central value for the lapse function N{0) between General 
Relativity and Brans-Dicke Gravity. Regardless of the node-number or value for A, there are essentially no differences 
between the two theories. If we take even larger values for A then the differences are even smaller. However, if we 
take smaller values for lobd then the differences get bigger (see Ref. for complete details on the n = case 
for ujbd — 6). This increase is natural since it is well-known that smaller values oi ujbd generally result in larger 
differences with General Relativity. It is thus not too surprising that when ujb d = 600 we cannot generate significant 
Scalar- Tensor effects. 

On the other hand, the quadratic coupling depends explicitly on the gravitational scalar field. Hence, even if the 
coupling is made to satisfy observational constraints at spatial infinity, it is still possible that the gravitational scalar 
can be made large enough inside a boson star that significant deviations from General Relativity can be produced. 
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We have constructed configurations for the quadratic coupUng case with k = 0.38, which is very close to the limiting 
value of K = 3/8 for this theory. Again, we take <i>c = 0.15. The quantitative results for <I>, N and A are virtually 
indistinguishable from the previous Brans-Dicke case, and thus General Relativity. The plot of the gravitational scalar 
ip, however, is different (see Figs. 4a, 4b, and 4c) in that it is positive, but still maintains the same basic shape as the 
Brans-Dicke case. The masses, particle numbers, radii, and central value for the lapse function A'^(O) for zero-, one-. 



and two-node solutions for $c = 0.15 and A = 0, 10, 100 are listed in Table III. Comparison with the previous two 



tables shows no important differences in any of the values. Also, if we take larger values of n the differences are even 
smaller than those for k = 0.38. It is because the observational constraints force the gravitational scalar to be small 
at spatial infinity that this coupling is unable to generate any significant differences with General Relativity. That is 
we see in Figs. 4a-c that the field equations are able to change the gravitational scalar by factors of order unity, but 
not factors of, say, a hundred. 

The second way in which we are to extract differences between General Relativity and Scalar- Tensor Gravity is to 
consider sequences of zero- node configurations, where individual members in a sequence are delimited by their value 
for $c- In Fig. 5 we have a plot of M, Np vs $c, for A = 0, 10, 100. For each value of A we see that there is an absolute 
maximum mass, and number of particles, that any boson star can have, exactly as in General Relativity. (Similar 
results for the mass were obtained by Gunderson and Jensen [T^ for the Brans-Dicke coupling with ujbd = 6, as well 
as the couplings considered by Torres fl^ .) Not unexpected, the value for the maximum mass and particle number 
are not changed significantly from the General Relativity results. In Fig. 6 we have plots of the binding energy versus 
the particle number. We note that all branches starting from <I>c = have negative binding energies, just like General 
Relativity. The corresponding plots for the quadratic coupling (for p = 10, and k = 0.38) contain essentially the same 
features as Brans-Dicke Gravity. 

So, we see that the observational constraints are so restrictive that no significant deviations from General Relativity 
can be produced. However, we can at least extract some generic behaviour pertinent to the gravitational scalar, and 
for boundary conditions consistent with observational constraints, from Figs. 7 and 8. They are the Brans-Dicke and 
quadratic coupling graphs, respectively, of M vs ip{0) for a sequence of zero-node configurations having A = 0, 10, 100. 
Each curve is parameterized by $c, and starts at M = 0. As $c gets bigger, then i^(0) decreases and M increases. 
At precisely the same value of $c where the maximum mass occurs, M reaches its maximum value in Figs. 7 and 8; 
also, the local minimum M in both Figs. 7 and 8 corresponds to the first local minimum in Fig. 5. 

The first point to be grasped from these plots is that the central value of the gravitational scalar has an absolute 
minimum value. The second is that there is a region of the graph where the same value for (p{0) corresponds to two 
different mass values. There is even the "crossover" point where two configurations have exactly the same mass. (It 
can be inferred from the discussion below that one configuration is stable, and the other is not.) The final point is 
that there is the asymptotic region where the mass M oscillates slightly in value, but (p{0) continues to grow. This 
is certainly surprising, since one would expect a growing (p{0) (which is related to the Newton gravitational coupling 
between two particles) to force the mass to change monotonically. 



V. STABILITY ANALYSIS VIA CATASTROPHE THEORY 



Catastrophe theory is a relatively new mathematical tool to explain a variety of changes of state in physical systems 
[^|. Rigorous theorems have been established demonstrating its validity. It is particularly adept at extracting 
discontinuous changes of state when gradual changes occur in the system parameters. It has been used in General 
Relativity to analyze the stability of black holes non-abelian black holes |2^, and boson stars Most 

importantly for the present work, is that application of catastrophe theory to boson stars [ p^ delivers the same 
results for changes in stability that are obtained from dynamical, numerical analysis [ p3| . 

Catastrophe theory can be used to determine stability if the system under analysis develops bifurcations, or cusps, 
when curves existing in its so-called equilbrium space are projected into its lower-dimensional control parameter space. 
The variables for the equilibrium space are given by a potential function, control parameters, and state variables. The 
variables for the control parameter space are the potential function and the control parameters. For the problem at 
hand, the equilibrium space is three-dimensional and the control parameter space two-dimensional, with the binding 
energy M — mNp taken to be the potential function, the total particle number Np chosen to be the control parameter, 
and the central value of the matter scalar field taken to be the state variable. An equivalent variation on this 
choice for the equilibrium space and the control parameter space is to replace the total particle number with the mass. 

The appearance of a cusp in a curve's projection is not enough to determine a change in stability. Rather, if it 
is known that a system is stable along one branch of a curve in the control parameter space, then it will become 
unstable on the next branch formed at a cusp. The trick is in determining stability along the previous branch. Like 
the case of General Relativity 11^, we can assume that boson stars along the very first branches in Fig. 6 (i.e., those 
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that start at Np = and end at the cusp labled Bi) are stable. They are the only ones composed entirely of negative 
binding energy states, they contain the flat-space limit, and furthermore, results of a dynamical, numerical analysis 
demonstrate that they are stable against small perturbations p^ . 

Fig. 9 is the typical equilibrium space for the Brans-Dicke coupling. The projection into the (7Vp,<i>c) plane at 
M — mNp — gives the same curve for the particle number as in Fig. 5. The projection into the {Np, M — mNp) 
plane at the maximum <i>c value gives the same curve as in Fig. 6. It is this second projection that yields the control 
parameter space. Hence, the cusp labeled Bi in Fig. 6 represents a change in stability, since the configurations on 
the branch leading up to Bi are all stable. At the cusp we have $c = 0.09 and a mass of M = 2.255/G*to, which 
is not significantly different from the corresponding General Relativity result (see The cusp at B2 does not 

represent a change in stability since the binding energies near this point are positive; likewise, for the cusp at i?3. 
As one might expect, the maximum stable masses for A = and 10 are not significantly different from their General 
Relativity values. Neither are the quadratic coupling results, since the equilibrium space and control parameter space 
are virtually identical to what is depicted in Fig. 9. 



VI. CONFIGURATIONS IN THE EARLIER UNIVERSE 



We have seen that the attractor nature to General Relativity of the quadratic model extends also to a sequence of 
equilibrium boson star configurations. Therefore, we repeat here the same calculations as the previous section but for 
earlier eras for the universe according to the quadratic coupling model. Goin g to earlier eras means taking smaller 
values for p (which, recall, is related to the redshift via z ~ e^°~P — 1) in Eq. ( 3.17 ). Our interest will be focused on 
the maximum allowed stable mass. 

A not so trivial point is that all our calculations must be made within the same Scalar- Tensor theory. What this 
means is that we are allowed to choose the additive constant to the coupling function a{Lp) only once, and not change 
it when we change the value for p. We will stick with the choice that makes a((poo) = at p = 10. Thus, a(<poo) is 
not zero for any other value o f p. T his complicates the determination of the ADM mass, sin ce fo r p 7^ 10 it can no 
longer be obtained from Eq. (3.37). Fortunate ly, w e can use t he energy functional of Eq. (3.10) to determine the 
total mass, once the rescaling relations of Eqs. ( |3.28 ) and ( 3.39 ) are properly taken into account. 

In Fig. 10 we have a plot of M vs We only go far enough so that the first extremum can be clearly identified, 
which, as we know from the previous section, corresponds to the maximum stable mass. Clearly, the main effect is to 
decrease the mass, as suggested by Torres |l^. If we go to early enough cosmological eras {p < 5), however, we find 
that the binding energy is positive, for all values of $c (see Fig. 11). Furthermore, the cusp that appears at Bi in 
Fig. 6 is clearly gone for p < 5. This indicates, then, that no stable Scalar- Tensor Gravity boson star can be formed 
during early enough cosmological times. 



VII. CONCLUDING REMARKS 



The main point of this work is to see if significant Scalar- Tensor Gravity effects can be generated in equilibrium 
state boson stars, even when the fairly restrictive observational constraints are imposed. In contrast with earlier works 
[^,0, we constructed not only zero-node solutions, but also one- and two-node solutions. We also illustrated some 
of the interesting features of the gravitational scalar for the zero-node case, and used catastrophe theory to determine 
the maximum stable mass for the zero-node states. Finally, we confirmed Torres' |Q suggestion that masses decrease 
below corresponding General Relativity masses for boson stars formed in earlier eras of the history of the universe. 

However, our results also show that stars formed too early are unstable. This suggests, for instance, that a 
scenario whereby boson stars are formed in some cosmological phase transition will not work in Scalar- Tensor Gravity. 
Moreover, there is every likelihood that any type of stellar object will suffer from the same instability. For instance, 
Harada has recently shown that there are a range of values for the coupling for which perfect fluid stars are 
unstable, although no explicit connection with the ambient cosmological conditions is made. 

We did not apply catastrophe theory to the higher-node solutions, because we know, from numerical evolutions 
p8[ , that they are all unstable in General Relativity. Dynamical 1-D evolutions are being performed to see if the 
higher-node solutions are generically unstable in Scalar- Tensor Gravity. As mentioned earlier, the code developed 
for these evolutions confirms that zero-node configurations on the initial branch of Fig. 6 are stable. It is also 
being investigated whether the zero-node states on the other branches collapse, or disperse. So far, the results do 
substantiate the conclusions obtained here using catastrophe theory. All these results are currently under preparation 
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TABLE I. Listed are some sample values for General Relativity, for ADM mass (modulo 1/G»m), particle number (modulo 
1/G«m^), radius (modulo 1/m), central value for N/(rn/Q,), and m/Q, for $c = 0.15. 
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TABLE 11. Listed arc some 


sample values 


for Brans-Dickc, for ADM mass (modulo 1/G, 


m), particle number (modulo 


1/Gtm^), radius (modulo 1/m) 


, central value for N/(rn/Q), and 


m/Q, for u)BD = 600 and $c = 
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TABLE in. Listed are some sample values for the quadratic coupling, for ADM mass (modulo 1/G*m), particle number 
(modulo 1/G»m^), radius (modulo 1/m), central value for N/{m/Q), and m/Q for k = 0.38 and 3>c = 0.15. 
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Figure captions 



Sample configurations of equilibrium state boson stars in the Brans-Dicke theory. Matter scalar field $ is plotted in 
the case of node=l (a), 2 (b), and 3(c), respectively. Solid, dotted and three-dot-dash lines are for A = 0, 10 and 100 
case, respectively. 

Fig§ 

Lapse N and the metric A for the solutions plotted in FigjlJ. The lapse N has been re-scaled to its original form so 
that its asymptotic value is 1. 

FigjI 

Gravitational scalar field (p for the solutions plotted in Fig.|l|. 
Fig^ 

Gravitational scalar field (p for the quadratic coupling model. The value for $c is the same as in Fig.yJ. 
Fig§ 

Mass and particle numbers versus central matter scalar field in the Brans-Dicke theory. 
Fig^ 

Binding energy M — mNp versus particle number Np in the quadratic coupling model. 
FigJ3 

A sequence of equilibrium solutions in the Brans-Dicke theory. Mass versus central gravitational scalar field value is 
plotted. 

Fig§ 

A sequence of equilibrium solutions in the quadratic coupling model. Mass versus central gravitational scalar field 
value is plotted. 

FigQ 

Equilibrium configurations in the equilibrium space. Binding energy M — mNp, particle number Np and the center 
matter scalar field $c are taken as the potential function, control parameter, and state parameter, respectively. 
One equilibrium sequence line and three projected lines of it onto 2-parameter planes are shown. The plot on the 
{Np, M — mNp) plane is identical with Fig.^. 
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Fig 

Changing the boundary condition for Lp^o via the cosmological solution of the quadratic coupling. The sequences of 
equilibrium solutions for cosmological time parameter p = 10(solid line), 6 (three-dots-dash line) and 5(dotted line) 
are plotted. The solid lines are identical with those in FigjH 



Fig.ll 



Binding energy vs. particle number Np for different cosmological time p — 10, 9, • • • , 1. The case for p = 10 is the 
same as Fig.|. Only the range < ^ 0.5 is plotted. 
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FIG. 1. Sample configurations of equilibrium state boson stars in the Brans-Dicke theory. Matter scalar field $ is plotted 
in the case of node=l (a), 2 (b), and 3(c), respectively. Solid, dotted and three-dot-dash lines are for A = 0, 10 and 100 case, 
respectively. 
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FIG. 2. Lapse A'^ and the metric A for the solutions plotted in Fig.|l|. The lapse A'^ has been re-scaled to its original form 
so that its asymptotic value is 1. 
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FIG. 3. Gravitational scalar field for the solutions plotted in Fig.^. 
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FIG. 8. A sequence of equilibrium solutions in the quadratic coupling model. Mass versus central gravitational scalar field 
value is plotted. 
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FIG. 9. Equilibrium configurations in the equilibrium space. Binding energy M — mNp, particle number A^p and the 
center matter scalar field "I>c are taken as the potential function, control parameter, and state parameter, respectively. One 
equilibrium sequence line and three projected lines of it onto 2-parameter planes are shown. The plot on the {Np,M — mNp) 
plane is identical with Fig.^. 
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FIG. 10. Changing the boundary condition for ifao via the cosmological solution of the quadratic coupling. The sequences of 
equilibrium solutions for cosmological time parameter p = 10(solid line), 6 (three-dots-dash line) and 5(dotted line) are plotted. 
The solid lines are identical with those in Fig.B 
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